Generating Primordial Black Holes Via Hilltop- Type Inflation Models 



0^ 
O 
O 
(N 

D 

in 



Laila Alabidi^'B and Kazunori Kohri^'Q 

'Astronomy Unit, School of Mathematical Sciences, Queen Mary University of London, 

Mile End Road, London, El 4NS, United Kingdom 
'^Department of Physics, University of Lancaster, Lancaster, LAI 4YB, United Kingdom 

It has been shown that black holes would have formed in the early Universe if, on any given scale, 
the spectral amplitude of the Cosmic Microwave Background (CMB) exceeds Vq ^ 10"'*. This value 
is within the bounds allowed by astrophysical phenomena for the small scale spectrum of the CMB, 
corresponding to scales which exit the horizon at the end of slow-roll inflation. Previous work by 
Kohri et. al. (2007) showed that for black holes to form from a single field model of inflation, the 
slope of the potential at the end of inflation must be flatter than it was at horizon exit. In this work 
we show that a phenomenological Hilltop model of inflation, satislying the Kohri et. al. criteria, 
could lead to the production of black holes, if the power of the inflaton self-interaction is less than 
or equal to 3, with a reasonable number or e— folds. We extend our analysis to the running mass 
model, and confirm that this model results in the production of black holes, and by using the latest 
WMAP year 5 bounds on the running of the spectral index, and the black hole constraint we update 
the results of Leach et. al. (2000) excluding more of parameter space. 
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I. INTRODUCTION 



The temperature anisotropies of the Cosmic Microwave 
Background (CMB) have been measured by WMAP on 
angular scales down to 6 ^ 0.3°, whereas they have yet to 
be measured to an effectual degree of accuracy on scales 
9 < 0.3° In fact, CMB data naively allows for a 
very large spectrum on these scales, i.e. a spectrum a 
few orders of magnitude above ~ 10~^, evaluated 
at horizon exit. One can, however, place an upper limit 
on the smaller scale spectrum by taking into account as- 
trophysical and cosmological constraints on black holes 
0, II, E H H S B B [13, Tni • if this spectrum is in fact 
close to this upper limit, then the situation allows for 
'large' fluctuations; large enough to collapse into black 
holes, known as Primordial Black Holes (PBHs) (c.f. 
[H, [il, M, [H, [H, [O, [3). Since the uncertainties in 
the primordial spectrum at such small scales are dom- 
inated by instrumental noise [l|, as opposed to cosmic 
variance which is the dominant source of uncertainties at 
larger angular scales, it may be that these uncertainties 
can be reduced in future surveys. Therefore, the question 
of whether the spectrum of perturbations on small scales 
is large enough to form PBHs is one that can, in theory, 
be answered. 

The 9 > 0.3° spectrum has been used extensively in 
discrimina ting between models of inflation (c.f. [l^, [20l . 
HH, [H, m, Hj, [2^). These analyses are based on the as- 
sumption that the anisotropies in the CMB, and hence 
the origin of large scale structure, are sourced by quan- 
tum fluctuations in a scalar field during or straight after 
inflation. As such, we assume that the fluctuations which 



sourced the PBHs are also generated during inflation, 
specifically towards the end of inflation, as it is during 
this epoch that the small scale fluctuations exit the hori- 
zon. Measuring the 9 < 0.3° spectrum will therefore not 
only probe generic signatures of inflation 0, [1^ but also 
act as an indicator for the shape of the inflationary po- 
tential [l^, [131 • In this paper we aim to exploit the latter 
purpose of measuring the small scale spectrum, and in- 
vestigate whether single fleld models of inflation can lead 
to the formation of PBHs. ^ 



tion 



To be specific, the primordial black hole (PBH) condi- 
[i, [l3, [H, \M, [HI can be expressed as: 
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where the subscripts e and * refer to the end of infia- 
tion and horizon exit respectively. Assuming a constant 
spectral index then: 
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where we used (jl.ip in the last semi-equality. AA^ ~ N 
refers to the number of e— folds that elapse from the time 
when scales of cosmological interest leave the horizon till 
the end of infiation. 

Taking a standard value of A^ ~ 60, it then follows 
from (|I.ip that ~ 1.3. This is beyond the upper limit 
of 0.993 allowed by the recent WMAP data, at 95% con- 
fidence limit with no tensor modes or running. Therefore 
we consider the variation in the spectrum up to second 
order, and assume that the spectral index depends on 
scale, i.e. V oc /^"sC^)-!^ then: 
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The Number of e— folds 
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where we used Ink = ln{aH) ~ N in the second step. 
Taking Ug ~ 0.95 and iV = 60 then requires n'^ ~ 0.0061 
to satisfy the WMAP bounds and produce Primordial 
Black Holes. 

We now wish to rewrite the Primordial Black Hole 
Condition in terms of the slow roll parameters. 

Recalling that the spectrum can be written in terms 
of the potential V and the slow roll parameter e — 
mli{V'/Vf/2, then [H: 
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Defining a new parameter B = Ce/e*, and combining 
equations (jl.ip and (|1.5p we find that the condition for 
PBH formation, without violating the aforementioned as- 
trophysical and cosmological bounds is: 



6-10" 



(1.6) 



An 'absolute' upper bound on the spectrum is given by 
[lo| as V^^^ 10~^, which translates to a lower bound 
B > 10-^' 

Equation (|1.6p tells us that for an inflationary potential 
to lead to the production of PBHs, its' slope must flatten 
towards the end of inflation. This shape is satisfied by a 



phenomenological model akin to the one analysed in 43 1 , 



and also the running mass model, first introduced in 
In these types of models, we require that the inflaton ini- 
tially be sitting at around the top of the hill, that is near 
a local maxima. This condition can be considered natu- 
ral [i^l from the viewpoint of eternal inflation and 
can be understood as follows: via some mechanism, be it 
quantum tunnelling or an inhomegenous pre-inflationary 
universe, the inflaton will somewhere, at some time find 
itself sitting at the top of the potential, at which point 
the universe will start to inflate. As long as the inflaton 
is undisturbed, the universe will inflate indefinatley, and 
can end up volumetrically dominating the universe. Since 
this process can lead to an indefinatley large volume, then 
even if there was the smallest probability that inflation 
were to start, it would [i^, |43]. Within this patch, quan- 
tum fluctuations in some sub-patches displace the infla- 
ton from its vestige causing it to roll either to the left 
or the right, and ending inflation in those reg ions, while 
overall, the patch continues to inflate [i^, We intro- 
duce these models in sections (|TT| and (jHip respectively. 
We present our results in section ((TV| and discuss them 
in section (fV)) . 



In this paper we use the duration of inflation, otherwise 
known as the number of e— folds iV, as a discriminator. 
It is defined as the ratio of the scale factor a at the end 
of inflation to a at the 'beginning' of inflation: 



iV = In 
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where the final semicquality comes from the slow roll 
approximation. 

To get a proper handle on how long inflation lasted 
from the time of horizon exit, one needs a complete his- 
tory of the Universe. At present though, we do not have 
an agreed upon mechanism of reheating. Therefore, one 
assumes an instant transition from inflation to a radia- 
tion dominated universe, and gets the bounds [49l |: 



10 < N < 110 



(1.8) 



The lower bound comes from the assumption that Nucle- 
osynthesis is well bounded, and the upper bound assumes 
that the universe underwent a few bouts of 'fast' roll in- 
flation. We do note that these are extreme bounds, and 
that the limits = 54 ± 7 are more widely acceptable 
(c.f. SHSi). 



B. Slow Roll and Cosmological Parameters 

The lowest order slow order parameters are given by: 
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which we then use to compute the spectral index and the 
running: 
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We use the bounds on cosmological parameters given 
by the WMAP year 5, Baryon Acoustic Oscillations and 
Supcrnovae data sets [sot : 



0.939 < < 1.109 
-0.0728 < n' < .0087 



(1.15) 
(1.16) 



where a zero tensor mode was assumed in the prior. This 
prior is reasonable since we are considering small field 
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models, characterised by a field variation that is smaller 
than the Planck mass A</) < mp\. In these models the 
gravitational wave signature will be small (sij . and by 
small we mean well below the sensitivity of WMAP5 pa- 
rameter estimation, so we do not calculate the associated 
parameter. 



II. THE TREE LEVEL POTENTIAL 



We consider the potential of the form: 
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where p > 2 and q > p, plotted in Fig. [TJ The case 
p — 2 and q > 4 can be generated from a flat di- 
rection in the Minimal Supersymmetric Model (MSSM) 
[13, m, [13, [H, [li, [13. In this case q is the non- 
renormalisable operator that depends on the flat direc- 
tions. One can also motivate the parameter range p = 2 
and q > 3 in [13, [H, [1§| , in this scenario the inflaton 
higher order terms are not Planck suppressed, and one 
gets a lower energy scale inflation, on the order of the 
TeV - GUT scale. We also maintain (p < mpi, a realis- 
tic bound from an effective particle physics perspective, 
which demands that one not consider mass scales larger 
than the largest naturally occurring scale, in this case 
the Planck mass. Then due to the Lyth bound, the grav- 
itational wave contribution of this model is negligible, 
regardless how long inflation lasts. 

In this setup we require that the potential at the end of 
inflation be flatter than it was at the time of horizon exit, 
so the inflaton must roll towards the origin. We denote 
the inflaton value at the maximum of the potential as 
(j)m, at horizon exit as 0* and at the end of inflation as 
(j)e. We impose the conditions: 



theoretically motivated setup, the running mass model 
[43,[60,[6l|,[63,[6l,[63,[6i,[6i,|63 which has the potential: 
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(3.1) 



where the mass of the inflaton is scale dependent and can 
be expressed as: 
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where is the mass of the inflaton squared, Aq is the 
gaugino mass squared in units of JTipi, and a is related 
to the gauge coupling. 

The potential can then be written as: 



V = V„ 



1 9 



mpi / 2(1 + Q;ln(0/mpi))^ \mp\ 

(3.3) 

and Bq ~ IJ.Q + Aq, with inflation occurring in the regime 
4> ^ mp\. This potential has the shape in Fig.[2l and the 
parameters have theoretically motivated constraints [6l[ : 

1 < /io < 0(10) 
< ^0 < C(10) 

10"^ <a< 10"^ (3.4) 



A. Linear Approximation 

The linear approximation of the running mass poten- 
tial is given by [23, [ll] : 



^-1- Y(M*+cln(0/(/)„O) 



(3.5) 



0* < (Pm 
<(f>e <(j>o 

where (j)o is the inflection point, 
given by (|2.ip : 



(2.2) 
(2.3) 

and (j)o are then 
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which is basically equation (j3.3p expanded about the 
maximum of the potential. The terms c and fi^ are re- 
lated to the theoretical parameters Aq, /Iq and a by: 
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Since deflnes the maximum of the potential then /ij = 
— c/2, this allows us to write: 



III. THE RUNNING MASS MODEL 

With the exception of p = 2 and integral values 3 < 
g < 9, the previous model is a phenomenological one. 
However the shape of potential does appear in a more 
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FIG. 1: For the figures on the left, increasing the variant decreases A'', and the opposite is true for the variants defining the 
figures on the right. 




FIG. 2: Plot of the running mass model. The solid line rep- 
resents a larger a than the dashed-dotted line. 



where a = — cln(0e/</'m). 

This approximation is only valid near the maximum. 
We neither expect nor get reasonable values of <^e using 
this estimate. However given 0* and </>£, we found that 
the linear approximation for N in p.7p appears consis- 
tent with the numerical calculation. 



IV. RESULTS 
A. Hilltop 

For the tree level potential, Fig. [3] depicts the depen- 
dence of N on "p and g, we note that N increases sharply 
as p increases, as expected. We then searched for the 
range of p and q parameters that satisfy the bounds 
10 < < 110 and found that the condition for PBH 
formation within this range is 2 < p < 3 and p < g < 4. 
Larger values of p lead to A'' ^ 110 while larger values of 
g/p do not satisfy the WMAP bounds, since they steepen 
the potential and result in an increased n' . We also found 
that p ~ 2 and g < 3 (note <? 7^ 3) places ]>S in the range 
TV = 54± 7. 

We then plot the dependence of iV on r\p and 77^ 
in Fig. [4l As expected N increases for decreasing r\p 
and vice versa for r\q. Once we have filtered out the 
reasonable values of A'^, we find that within the range 
\J]p^'^q\ = {0,1} PBH formation will occur. It seems 
that the stronger constraints for PBH formation with 
10 < TV < 110 come from p and q. 

Next we consider the case of defining by the condi- 
tion that n(«i,) = 0.95 and by iV = 60 and A^ = 100. 
Fig. [5] shows that for both A^ = 60 and iV = 100, the pa- 
rameter ranges satisfy the WMAP bounds and the PBH 
constraints. However, for A'^ = 60 integral values of p 
and q do not lead to the formation of PBHs, while for 
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FIG. 3: A plot of the maximum and minimum values of 
log(log(A'^)) versus p and q, from a range values of r/p and 




FIG. 4: A plot of the minimum values of log(A'^) versus r)p 
and r)q, from a range values of p and q. 

N — 100, the parameter set {p, q} ~ {2, 3} does. 

B. The Running Mass Model 

In this analysis we considered the parameter ranges 
in p.4p . for which we find that this model generates a 
large range of N values. In Fig. [6] we plot the allowed 
parameter space for 10 < < 110 and iV = 54 ± 7. 

From Fig.[7]wc find that in order to avoid the overpro- 
duction of primordial black holes, n > 1 would be ruled 
out for iV = 45, and this bound is strengthened to rul- 
ing out n > 0.95 for N = 60. These bounds are slightly 
stronger than those found by [67j who rule out n > 1.1 
for N = 45. 

However, there is a discrepancy between our spectral 
index contour lines and (67| , which we found was resolved 



by evolving our system an extra ~ 7 e— folds. Via a pro- 
cess of elimination, we think this may be due to the fact 
that Ref. [HI] solved the background equations numer- 
ically without resorting to the slow roll approximation, 
while using the extended slow roll formalism of [1^ to 
solve for the perturbations. 

Either way, as our method underestimates the allowed 
parameter range for each n, then using ref. [67| method 
would strengthen our conclusions i.e. our bounds are con- 
servative. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we utilised the spectrum on scales 9 < 
0.3°, corresponding to the end of inflation, to further the 
field of inflation model discrimination. The spectrum on 
these scales has yet to measured, but future CMB surveys 
such as the PLANCK mission may constrain its' value. 
Astrophysical phenomena determines the upper bound 
to be ~ lO"**, corresponding to the criteria for the 
formation of primordial black holes (PBHs). In terms of 
the e slow roll parameter, this means that the value of e at 
the end of inflation must be much smaller than it's value 
at horizon exit p.6|l . The models of inflation that satisfy 
this condition must therefore exhibit the unique property 
of having a flatter slope towards the end of inflation. So 
far, we only know of generic Hilltop models of inflation 
that fulfil this criteria. 

We have investigated whether these generic models of 
Hilltop inflation would lead to the production of primor- 
dial black holes with a spectrum Vq^ ^ 10~^. We found 
that within the range of parameters allowed by the lat- 
est WMAP data, the Hilltop model (HHD would lead to 
the formation of PBHs without violating astrophysical 
bounds for p < 2.5 and g < 3 if > 60, and for p ^ 2, 
2<g<c3ifiV> 40. Integral values of p and q, which 
have some theoretical motivation, only lead to PBH for- 
mation within the bound (|1.6p for p = 2 and 9 = 3, with 
60 <C iV < 100. In all cases it seems that near maxi- 
mal running is required. If, however we were to allow 
^ 110 the range of parameters that would lead to 
PBH formation would be extended. 

The allowed parameter range for the production of pri- 
mordial black holes with Vq^ ^ 10~^ in the running mass 
model is again dependent on N as can be seen from Fig.[Sl 
We find that for a = 0.01 and /Xg > 1.1, black holes could 
form after iV < 47 e— folds, and therefore before what 
can be considered a 'reasonable' end to inflation. This is 
problematic on two counts, if we assume that the PBHs 
formed prior to the end of inflation, then this could lead 
to the overclosure of the universe (c.f. [6§|). On the same 
note, we know that on CMB scales the spectrum is too 
small to support PBH production. On the second count, 
assuming that the formation of the PBHs coincided with 
the end of inflation, then the arguments we presented in 
section p Ap apply. Thus, using as a discriminator we 
rule out ^0 > 3 and /Iq > 1.1 for a = 0.01, Aq > 6 and 
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^0 > 2.4 for a = 0.005, and > 5 and > 8.75 for 
a = 0.001. As we mentioned in the text a = 0.1 is ruled 
out on WMAP consistency grounds. 

This model has also been analysed by [T^I, in which 
the authors use neutrino and 7— ray background data 
to constrain the PBH mass spectrum, which determines 
the spectral index on small scales /c ~ 15 Mpc~^. Then 
assuming that the running mass model is correct they 
reconstruct the power spectrum, finding that on these 
small scales the spectrum is highly sensitive to the run- 
ning of the spectral index. Combining these two pieces 
of information they get bounds on and n'^ , which turn 
out to be inclusive of the WMAP limits on these pa- 
rameters. That is, by using a different approach to ours 
(70I I conclude that the running mass model is consistent 
with WMAP while avoiding PBH over-production, con- 
gruously with our findings. 

Finally, we note that our generic results are consis- 
tent with the findings of [7lj. They tackle the question 



of PBH production utilising the arising constraints to 
derive bounds on the cosmological parameters, and con- 
clude that the PBH constraint is 'strongly' dependent on 
N and the spectrum at the end of inflation. Character- 
istics exhibited by our specific models. 
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FIG. 5: Plot of log(S) versus log(n') for the Hilltop model with TV = 60 (figure on left), A'' — 100 (figure on right) and Us = 0.95. 
The hatched region is excluded, representing log(n') > —2 and log(Z?) < —8. The region log(B) > —6 does not lead to the 
formation of PBHs, and is represented by the tan colour in the figure. PBHs can form in the region —8 < \og{B) < —6 without 
violating astrophysical or cosmological bounds, and is represented by the light orange region. The yellow dots correspond to 
{Pi q} = {3, 4}, the green dots to {p, q} — {2, 3} and the blue dots to {p, q} = {2, 2.5}. 




FIG. 6: Contour plots of the number of e— folds produced in the running mass model for three values of the gauge coupling 
a = [0.001,0.005,0.01]. We found that a = 0.1 did not satisfy the WMAP bounds on Hs and n . We have fihered out the 
allowed parameter space for 10 < N < 110, and coloured it in shades of green. The dashed regions in each plot correspond to 
the more 'reasonable' bound A'' = 54 ± 7 
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FIG. 7: In these plots we fix and a, plotting contour lines of the spectral index (dashed) and log(S) (solid). Note that the 
contour lines do not exactly match Ref. [g^l, an anomaly that we discuss in the text. Parameter space below log(S) ~ —8 is 
excluded. 



